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We consider a stochastic model of investment on an asset of a stock market for a prudent 
investor. She decides to buy permanent goods with a fraction a of the maximum amount 
of money owned in her life in order that her economic level never decreases. The optimal 
strategy is obtained by maximizing the exponential growth rate for a fixed a. We 
derive analytical expressions for the typical exponential growth rate of the capital and 
its fluctuations by solving an one-dimensional random walk with drift. 



1. Introduction 

A large number of studies on finance have the main purpose to find the optimal 
strategy for a given kind of investment Sffl'!. These problems can be tackled 
by looking for simplified (but non-trivial) models which are able to describe the 
observed phenomenology and which can be, eventually, solved analytically. For an 
introductiouixLfinaiicial problems discussed from the point of view of the theoretical 
physics seeiwie. 

The optimal strategy is usually defined as the one which maximizes a given 
utility function which takes into account the risk. The use of utility functions 
introduce a high degree of arbitrariness since the particular choice depends on the 
subjective aversion to risk of the investor. This psychological arbitrariness can 
be removed if one considers that the rate of growth of the capital is an almost 
sure quantity on the long run. Therefore, the best strategy can only be the one 
which maximizes this rate i.e. the one which maximizes the expected logarithm of 
the capital. Any other strategy ends almost surely with an exponentially in time 
smaller capital BiH'B. 

The deep understanding of the reasons for the use of logarithmic optimization 
strategy comes from the pioneering Kelly's work 0. In his paper he proposes and 
solves a model where an investor uses a fraction I of her capital to buy shares of a 
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given asset at discrete time steps. It is assumed that the price of the shares can, 
at each time, double or vanish, so that the investor doubles or loses the fraction 
she has invested in it. It is assumed that the probability p of doubling is larger 
than 1/2, this is a reasonable assumption since the contrary situation implies that 
is better to keep the money in a free risk investment (the interest rate is supposed 
to be vanishing). 

If the investment is absolutely sure (p = 1), of course she will invest all the 
capital (I = 1) at each step. In this way after t steps her capital will be 2* times 
the original one. On the other hand, if the evolution of the share price is uncertain 
and she wants to maximize the expected value of her capital, she chooses the same 
strategy by investing a fraction 1 = 1. Obviously this is not the best one: in fact it 
is sufficient only one defeat to lose everything. 

Using arguments from the theory of probability, Kelly has shown that the correct 
quantity to maximize is the expected value of the growth rate x(t) of the capital: 
this quantity corresponds to the rate of transmission over a channel in information 
theory or to the Lyapunov exponent in dynamical systems and statistical mechanics 
of disordered systems. The value of x(t) f° r a particular sequence of investments 
fluctuates around the expected value (x(t)} and the fluctuations go to zero in the 
limit t — ► oo. 

If one introduces a random interest factor r t and/or a random gain factor Vt 
(in the original Kelly's work r t — 1 and v t = 2) the model is still trivial from a 
mathematical point of view and it is easy to find out the optimal strategy. This is 
due to the fact that the model can be written in terms of a multiplicative random 
process. For a discussion on optimal investment strategy of a multi-asset portfolio 
following Kelly's approach see 0. 

Recently Galluccio and Zhang i have considered a generalization of the previous 
model, where at each time step several kinds of investment are possible: a sure one 
(i.e. the bank) and the other ones risky (i.e. the stock market). They assume, as 
Kelly does, a simple behaviour for the market and find the values of the parameters 
which optimize the Lyapunov exponent. This model can be written as product 
of independent random matrices and it can he treated with standard perturbative 



methods 



■endent random matrices and it cap h 
i E3 or by constrained annealing Ejli-X 



In this paper we consider the case of an investment where there is a diversification 
between the stock market and permanent goods. A permanent good, such as a 
house, is characterized (at least as a first approximation) by a value that does not 
change in time, and it is not easy to convert in cash. The prudent investor wants, 
at least, to assure herself the same economic level for all her life (i.e. the capital 
invested in permanent goods cannot decrease) . In order to reach the goal she decides 
the permanent goods must equal a fraction a of the maximum total capital owned 
in the past. Therefore the capital cannot be less than this threshold and only the 
remaining part can be invested in the market. As a consequence the model describes 
the capital as a stochastic variable with memory. 

Let us briefly sum up the contents of the paper. In section 2, we introduce 



Optimal strategies for prudent investors 3 



our model as a modification of the Kelly's one and we give an interpretation of 
the parameters introduced. Section 3 is devoted to the analytical solution of the 
model. In section 4, we compare our results with the Kelly's ones. In section 5, we 
discuss some aspects of our model and its main features, in particular we consider 
the possibility of looking for time changing strategies of investment. 

2. The model 

A given asset in a stock market can be always modelled in absence of memory by 

W t+1 =F t (W t ). (2.1) 

where Wt is the capital at discrete time t and Ft is a random function, i.e. attach 
time t one chooses among different functions according to given probabilities liSEil 
The simplest case is F t = UtWt, where the u t are independent stochastic variables 
(e.g. they can assume only two values as in the case of Kelly). 

In the Kelly's model the investor keeps untouched a fraction l — l (with < / < 1) 
of its capital, while the rest is used to buy shares with two different results: or she 
doubles her investment, either she loses it. Therefore the capital at time t + 1 is 
given by 

W t+ i = (1 - l)W t + 1(1 + <T t )W t = (1 + la t )W t . (2.2) 
It follows that ut can be written as 

u t = 1 + lo t 
where the dichotomic random variable at 

{+1 with probability p 
— 1 with probability 1 — p 

describes the change of the share price. 

The growth rate of the capital at time t is %(i) = | In Wt. This quantity is 
random. Nevertheless for large t, because of the law of large numbers, it converges 
almost surely to the Lyapunov exponent 

A = lim X (t) . (2.3) 

t — >oo 

Let us notice that the optimal strategy proposed by Kelly consists in the maxi- 
mization of A (i.e. (In Wt)) and not of (Wt). Following the naive idea to maximize 
(Wt) one has I — 1 and (Wt) = (2p)' which is much larger (at large t) than the 
corresponding value obtained with the maximization of A. Nevertheless, the naive 
strategy is clearly wrong since at large t one has a probability close to 1 to lose ev- 
erything. On the contrary the maximization of A has a well established theoretical 
motivation in the law of large numbers. Basically one has that for almost all the 
realizations the quantity \(t) a t large t is close to its mean value A, i.e. the quantity 
x(t) is self- averaging. Then, since we are dealing with a multiplicative process, one 
has that the probability distribution of Wt is close to the log-normal one ExtLZI: 
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PW)* , 1 exp-^-^. (2.4) 

where 

A 2 = lim t{(- hxWt — A) 2 ) . (2.5) 



The (2^) holds for small values of (In Wt — Xt) /AviLon the contrary the tails 
depend on the details of the multiplicative process E30. Let us remark that the 
investor can have a small but finite probability to have a very low capital at time t. 

The optimization problem consists in choosing the fraction I of the capital that 
maximizes A given the probability p > 0.5; the result is l max = 2p — 1. 

Our aim is to modify the Kelly's model in order that the capital cannot become 
too small in any realization of the random sequence {at}. This time the investor 
decides to buy shares with only a part of the capital and she uses the other part to 
buy permanent goods. The value of the goods equals a fraction a of the maximum 
total capital she has owned in the past. 

The model can be written in the form: 

Wt+i = aW t + (1 + lcrt)(W t - aW t ) . (2.6) 

where 

Wt = max Wi. (2.7) 

{*<*} 

Let us stress that aWt is the part of the capital kept untouched (Wt+i is always 
larger than aWt), and I is the fraction of the available part Wt — aWt, risked at 
time t. In the following a will be considered a fixed parameter depending on the 
greediness (or the prudence) of the investor. The Kelly model is recovered for a = 0. 

Let us remark that the optimal strategy of the model ( |2.6| ) - ( |2.7D is, from a 
conceptual point of view, equivalent to the optimal strategy of the original Kelly's 
problem with a suitable utility function which takes into account the prudence of 
the investor. The model is then a non-markovian process for the single variable Wt. 
A typical realization of Wt and Wt is shown in Fig. 1. 

Moreover it is worth interesting that our model can be thought as a Markovian 
process of two variables (Wt and Wt) if we express the maximum capital owned in 



the investor's life (2/7) as: 



W t+ i = max{W t , Wt+i} = max{W t , (1 + la t )(W t - aW t ) + aW t } . (2.8) 



Equations (2.6) and (2.8) are in fact a random map of two variables of the form 

H '«w,m. (, 9) 



Wt+x ) L \w t 



Equation (2J3) can be thought as a product of random matrices of infinite order. 



This is clear if one consider the process ( |2.6| ) with where now Wt is 

W t = V t {1) = max Wi . (2.10) 

{i=t-i,t} 
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Fig. 1. A typical realization of the capital Wt (full line), of its maximum Wt (dashed line) and 
of the capital invested in permanent goods aWt (dotted line), for a = 0.7 and p = 0.75. The two 
arrows shows a process of type (3.12), where Wt is constant. 



or 



W t = V} j) = max 



{t-j<i<t} 



(2.11) 



It is easy to realize the model (2.6) with Wt given by ( 2.10 ) can be represented in 
term of a Markov process of order 1 (i.e. the state at time t depends from the states 



at t and t — 1). In a similar way using (2.11) one has a Markov process of order j 



As far as we know for a general problem like (2.9) there is not an Oseledec 
theorem El for the self-averaging of the quantity x(i). Nevertheless for our specific 
case we shall show in the next section that x(t) 1S self-averaging. 



3. The solution 



The process ( |2.6| ) and (2.7) can be thought as a sequence of independent processes. 
Each of them starts at time ij and ends at time ti+\ when Wt changes its value (i.e. 
Wt reaches a new maximum) so that, during this time, Wt is constant and equals 
the starting capital W ti (see Fig. 1). 

Between ti and t,_|_i the variable Wt — aW^ is multiplicative as in the Kelly's 



model, in fact the equation (2.6) reduces to 



W t+ i ~ aW u = (1 + lo t )(Wt - aW u ) . 



(3.12) 



where W ti plays the role of a constant. 
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We notice that the final value of the i process W* i+1 turns out to be propor- 
tional to its initial value Wt{- 



e 1 Wt, 



(3.13) 



where 7^ depends on all the {cr t } extracted during the time interval (ti,ti + i). 

This process ends when Wt i+1 becomes larger than W ti ( 7W > ) for the first 
time. Of course the time interval ATW = - f t is a random quantity, and it 
depends on the {ct} sequence. 

In this contest the global process Wt (2.6) can be expressed in terms of M 
independent Markovian processes r )^ 1 \..., r y^^ as: 

M 

W t = ^ J|e 7<,) ■ (3.14) 



where t equals the sum of the time duration of all the M Markovian processes 
t = Y^iLi N®. Let us remark that ( 3.14 ) is a product of independent random 
factors, this implies that %(£) reaches the value A at large t for almost all realizations. 

Let us stress that the i th process described by ( 3.12 ) is a one-dimensional random 
walk with positive drift (p > 1/2), in terms of the variable X)t'=t a t' ■ The process 
ends as soon as the random walk reaches an escape point that runs away with 
velocity ^py, where p is a monotonic function of I defined by 



ln(l - I) 
"ln(l + 



(3.15) 



with p > 1. 

Supposing that it happens with nS 1 ' defeats (or negative steps of the random 
walk), we find that 



In 



a+(l-a)(l-i) n (l + l) 1 



Let us notice that and are not independent but they must satisfy 



N {l) = 1 + n (l) + 



(3.16) 



(3.17) 



where the square bracket indicates the integer part. 

Let p(n) be the probability that the process 7 ends with n defeats; it can be 
written down as 

p(n) = C n (l - p) n p N -~ n . (3.18) 

where C n is the number of different ways to exit from the process with n defeats. 
The following recursive formula for C n holds (see Appendix) for n > 3: 



C n — 



N n -1 - 1 

n — 1 



n-2 



N n -i - N r 
n — r 



(3.19) 
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with initial values 



Co = Ci = 1 
C 2 = N t - 1 



(3.20) 



Equations (3.19) and (3.2C) represent a practical tool in order to numerically com- 
pute p(n). 

Let us notice that, fixed p, exists an interval of p (i.e. of I), such that the i-th 
process has a non vanishing probability to have an infinite time duration (it happens 
when the mean velocity 2p — 1 of the random walk is lower than the velocity 
of the escape point). This implies that the growth rate A of the capital is almost 
surely zero, since its evolution remains confined in a i-th process, with finite i. Of 
course it is a non optimal situation, and the following considerations are restricted 
to the more interesting cases with A > 0. 

In the limit M — > oo, because of the law of large numbers, the Lyapunov expo- 
nent A can be written as: 



A 



lim 



A/ 



M 



In 



(3.21) 



where the bar indicates the average according to the distribution p(n). 

Let us recall that the distribution of Wt is approximated by a log- normal djOh , 



and furthermore the variance (2.5) can be written as 



A = lim 

M— >oo 



(S,ti7 W -AS^i^ W ) 2 



(7n - xN n y 



(3.22) 



where the last result is obtained simply noticing that J2i=i is equal, for the 
law of large numbers, to M N + O(M^). 



4. Discussion of the results 

As in the Kelly's model we maximize the Lyapunov exponent with respect to I, to 
find the long time optimal strategy. We have computed A and its variance using 



equations (3.21) and (3.22). The probability p(n) is found out for any n smaller 

than an appropriated n so that J2n=o P( n ) — ^ — 10~ 8 , n is typically O(10 2 ). 

In Fig. 2 we report p max = p(l function of p for a = 0.6. We observe 

that p max is constant for some intervals of p. These plateaux correspond to rational 
values of p m ax- This feature of the model can be explained noting that the proba- 
bility p(n) is discontinuous around any rational p. As a consequence, the Lyapunov 
exponent A as function of p for a fixed value of p, has a cusp at any rational value 
of p. One of these cusps is a maximum of A in correspondence of the plateau of p 
(see, for instance, Fig. 3a). 

In this context, varying p one has only a rotation of the cusp, so that at different 
p corresponds the same value of p max that maximizes the Lyapunov exponent A, 
while for p out of the plateau one has a decreasing or an increasing cusp, like in Fig. 
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Fig. 2. pmax = p{lmax) (3.15) as function of p for a = 0.6. The plateaux correspond to rational 
values o{ p ma x- 



3b. The width of the plateau depends on a and it becomes larger when a increases 
and when p is integer. 

Let us restrict to the case of integer p. In order to simplify the notation, we use 
* to indicate a quantity computed at integer values of p. If we study the probability 



distribution of the defeats ( 3. IS ) for fixed values of a and p as a function of p we 
notice that, when p crosses an integer value p, the time durations {N n } ( 3.17 ) and 
the coefficients {C n } ( |3.19 ) change for every n. N n changes since it depends from 
the integer part of np, and C n since the number of paths of the random walk ending 
with n losses depends on the cases with less defeats. For the same reason both of 
them remain constant immediately on the right and on the left of p, at least till n 
is big enough to have negligible effects on the probability distribution p(n) . 

Following this idea we perform a linear expansion in I of A around its value on 
the cusp. Then we write 



X ± = A± + XfSl 



where A + and A are respectively the right and the left limit for / 
After some algebra one obtains 



A+ = Aq = A = (p - p(l - p)) \n(a + (1 - a)(l + I)) 



and 



1 - a 



p[2 + b{l + p)]-[l + l + pb]) 



I ± of A. 
(4.23) 
(4.24) 
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Fig. 3. Lyapunov exponent A as function of p for: 

a) p = 0.8 and a = 0.6, where the maximum is a cusp (p = 2); 

b) p = 0.6 and a = 0.6, where an increasing cusp (p = 1.5) is present but not in correspondence 
of the maximum. 
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1 - a 



al(l 



b[2p- (l + l)] 



where 



b = 



al(l - I) 



(4.25) 



(4.26) 



a + (1 - a)(l + 

The signs of A^ tell us when the cusp is a maximum of A as a function of I (i.e. 
when AT > and A^ < 0) . From (4.24) and (4.25) is easy to see that this happens 



when p is between p m i n and p max where 



Pmin — 



Pn 



1 + 1 + pb 
2 + 6(1 + p) 



l + l 



(4.27) 



(4.28) 



Some of the widths of these plateaux (i.e. p m ax — Pmin) are plotted in Fig. 4. 

It is interesting to observe that the plateaux disappear (p m ax — * Pmin) m the 
Kelly's limit (a — > 0) and p mQa; is the value for which the Kelly's model reaches 
its maximum when I — I. Then we notice from (4.23) that for this value of p the 
maximum Lyapunov exponent rescaled with the function 



rj(a,p) = ln(a + 2(1 — a)p) 



(4.29) 
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Fig. 5. Rescaled maximum Lyapunov exponent — (l m ax) (full line) as function of p for a = 0.6, 
compared with the Kelly's case (circles). 



is the same for every a. 

Finally we notice that the Lyapunov exponent computed at the optimal value 



I = lmax, rescaled with ( [4.29|) , is quantitatively independent of a, see Fig. 5. 

In addition the standard deviation A is proportional to the Lyapunov expo- 
nent when I = lmax- This is an exact result, at any a, when p = p max and it is 
qualitatively true for generic values of p, see Fig. 6. 

Let us remind that a good parameter for quantifying the strength of the fluctua- 
tion is the ratio R — A 2 /A. In fact, in the log-normal distribution both A 2 and 

a 2 

A have the dimension of the inverse of time, and therefore ^- is time independent. 
Since the A/ A is basically only function of p and not of a, one has 



A 2 /A\ 2 /A\ 2 T) rj 



R = —T- — I -T- I A ~ — I XKelly — RKelly < RKelly ■ 

A V A / Kelly V A / Kelly VKelly V Kelly 

(4.30) 

i.e. a reduction of the relative fluctuations. 
5. Conclusions 

In this paper we have considered a diversification of the portfolio between permanent 
goods and investments in a market. The model is a non trivial modification of the 
Kelly's one where only a part of the capital is allowed to be invested on the market. 
The investor keeps the remaining part as a security amount of money, equal to a 
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P 



Fig. 6. Relative deviations ^(l m ax) as function of p for a = 0.3. 

fraction a of the maximum capital owned in the past. In this way the investor 
avoids the possibility of losing almost all her capital due to a large fluctuation as 
it can happen in the Kelly's case: the parameter a can be thought as a measure of 
the investor's prudence. 

The small fluctuations of the capital around the typical value Wo e xt follow, 
at large t, a log-normal distribution and therefore they are well described by the 
typical exponential growth (or Lyapunov exponent) A and the deviation A from 
this quantity. We give explicit analytical expressions for both these quantities. In 
particular we obtain a decreasing of the relative strength of the small fluctuations 
(i.e. A 2 /A). 

An interesting feature of the model, from a mathematical point of view, is that 
the Lyapunov exponent is a continuous but not differentiablc function of the pa- 
rameters. This fact is particularly relevant when we look for the fraction I of the 
allowed capital (i.e. the capital that can be invested in the_market each time), which 
maximizes A. We observe a devil's stairs like behaviour 113 for p max as a function 
of the probability p. The existences of plateaux can be understood if one considers 
in more detail the A itself as a function of I at fixed p and a. The sizes of these 
plateaux can be computed, as an example we derive the width of the largest ones. 

We have found that the Lyapunov exponent, when is rescaled by a proper func- 
tion of the parameters, and the ratio between A and A, show a similar behaviour 
of the Kelly's case: the prudence's constraint we have introduced has basically the 
effect to rescale the exponential growth and the relative strength of the small flue- 
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tuations of the capital invested on the stock market according to equations (4.29) 
and ( p~30l ). 
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The study of non commutative multiplicative models of the stock market has 
the great advantage that they often can be analytically treated. In generic^ 
one can use powerful systematic methods to obtain good approximations 11311 
We belive they can be useful to understand problems close to the reality, such as the 
ones where the investor does not decide only once the best strategy to follow, but 
can change at each step his mind depending on the behaviour of the market. It can 
be shown that a fixed strategy is not the best one when we introduce correlations 
between successive times or different hedgings for prudent investors. 
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Appendix 

In this appendix we derive the recursive formula ( |3.19 ) for C n . The random walk 



(3.12) has a positive drift (p > 5) and it ends as soon as the capital exceeds its 
initial value. The total number of steps necessary for that is N n , where n is the 
number of negative steps (i.e. defeats). The stop is when n/N n is smaller than 
l/(p + 1) (with p > 1), so that N n = 1 + n + [np], where the square brackets mean 
integer part (see equation ( p. 17 )). 



The number of different ways C n can be computed noticing that the n negative 
steps have to appear before the (iV n _i)-th step, in order to avoid a premature 

interruption of the process with only n — 1 negative steps. This yields to ' 
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different combinations, but in this number are also included the cases of premature 
arrest with r negative steps in the first N r time steps, with < r < n — 2. Each 

of these cases leaves out an amount of combinations equal to [ n_1 r \ C r , 



where the combinatorial factor comes from the different ways that the remaining 
(n — r) negative steps have to appear in the interval [1 + N r , JV n _i]. 
It immediately follows the recursive formula ( 3.1 9| ) 



( " ■ In) ^ r=0 { n-r ) ° r " 
n-1 J ^=1 I n _ r I W 



which holds for n > 3, while for n = 2 one has 



C 2 = JVi - 1 
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:ond- 



The cases n < 1 can be trivially derived: 

C a = C\ = 1 
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